TO APPEAR IN THE ASTRO PHYSICAL JOURNAL 
Preprint typeset using LAT^X style cmulatcapj v. 5/2/11 
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ABSTRACT 

We present a composite light-curve model of the symbiotic nova PU Vul (Nova Vulpeculae 1979) that 
shows a long-lasted flat optical peak followed by a slow decline. Our model light-curve consists of three 
components of emission, i.e., an outbursting white dwarf (WD), its M-giant companion, and nebulae. 
The WD component dominates in the flat peak while the nebulae dominate after the photospheric 
temperature of the WD rises to logT (K) > 4.5, suggesting its WD origin. We analyze the 1980 and 
1994 eclipses to be total eclipses of the WD occulted by the pulsating M-giant companion with two 
sources of the nebular emission; one is an unocculted nebula of the M-giant 's cool- wind origin and 
the other is a partially occulted nebula associated to the WD. We confirmed our theoretical outburst 
model of PU Vul by new observational estimates, that spanned 32 yr, of the temperature and radius. 
Also our eclipse analysis confirmed that the WD photosphere decreased by two orders of magnitude 
between the 1980 and 1994 eclipses. We obtain the reddening E(B — V) ~ 0.3 and distance to 
PU Vul d ~ 4.7 kpc. We interpret the recent recovery of brightness in terms of eclipse of the hot 
nebula surrounding the WD, suggesting that hydrogen burning is still going on. To detect supersoft 
X-rays, we recommend X-ray observations around June 2014 when absorption by neutral hydrogen is 
minimum. 

Subject headings: binaries: symbiotic — nova, cataclysmic variables — stars: individual (PU Vul) — 
stars: late-type — ultraviolet: stars — white dwarfs 



1. INTRODUCTION 

Symbiotic novae are thermonuclear runaway phenom- 
ena occurring on white dwarfs (WDs) in binary systems 
that consist of a WD and a red giant (RG). Symbi- 
otic novae can be divided into two subgroups accord- 
ing to their spectral evolutions. The first group ex- 
hibit a long (several years) "supergiant phase," resem- 
bling an A-F supergiant when the star underwent an 
outburst. In the second group, a nebular phase begins 
almost immediately after the optical maximum, and a 
"supergiant phase," if there is, has a very short duration 
(jMurset fc Nussbaumer I fl994ft . The first subgroup in- 
clude AG Peg, RT Ser, RR Tel, and PU Vul. The second 
subgroup include V1016 Cyg, HBV 475, and HM Sge. 
It is, however, still unknown the reason why this differ- 
ence arises. Due to very long evolution-timescales (one to 
several tens of years or more) it is not easy to obtain ob- 
servational data of a whole period of the outbursts, such 
as dense and continuous photometric, spectroscopic, and 
multi-wavelength observations including UV and X-ray. 
Under these circumstances, it has been difficult to study 
symbiotic novae quantitatively compared with classical 
novae. 

Among symbiotic novae, PU Vul is a rare exception. It 
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is an eclip sing binary of the orbi t al period ~ 4900 days 
(13.4 yrl (iKolotilov et al. 1 119951: iNussbaumer fc Vogell 
[199^1 iGarhavich ||1996t IShugarov et al. 112011ft . During 
eclipses, different emission components are occulted dif- 
ferently. This offers a good chance for quantitative study. 
PU Vul outbursted in 1979 and we have dense opti- 
cal spectroscopic /photomet ric data as well as IUE/HST 
UV observations. Recently, iKato et al. I fj201 lh first pre- 
sented a theoretical model of PU Vul that reproduces 
the optical flat peak as well as the UV light curve, and 
estimated the WD mass (~ 0.6 Mq). Their model is, 
however, only for the light curve of the outbursting com- 
ponent (WD) and the other emission components were 
neglected. 

This paper presents a comprehensive model of emis- 
sion components of the WD, RG, and nebulae, based on 
new estimates of the temperature and radius of the hot 
component (WD), as well as the cool component (RG) 
derived from the two eclipses (1980 and 1994). Section[2] 
briefly introduces our evolution model of PU Vul. Using 
our theoretical light curves, we constrain the extinction 
and distance to PU Vul in Section [3J Section @] compares 
our theoretical light curves with our new observational 
estimates of temperature and radius of the WD. In Sec- 
tion [5j we analyze light curves of the two eclipses and 
obtain the binary parameters as well as the brightnesses 
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of the RG and nebulae. Using these values, we construct 
a composite light curve model of PU Vul in Section [6j 
Discussion and conclusions follow in Sections [7] and [8] 

2. MODEL OF WD COMPONENT 
2.1. Evolution of Nova Outbursts 

A nova is a thermonuclear runaway event on a WD. Af- 
ter the hydrogen shell flash sets in, the envelope on the 
WD expands to a giant size. After it reaches the optical 
peak, the envelope settles down into a steady-state. The 
optical magnitude decays as the envelope mass decreases 
while the photospheric temperature rises with time. The 
decay phase can b e followed by a quasi-static sequence 
(jKato et al. 1 120 111 ). We solved the equations of hydro- 
static balance, continuity, radiative diffusion, and conser- 
vation of energy, from the bottom of the hydrogen-rich 
envelope through the photosphere. The evolution is fol- 
lowed by a sequence of decreasing envelope mass. The 
time interval At between two successive solutions is cal- 
culated by At = AM env /(M nuc + M wind ), where AM env 
is the difference between the envelope masses of the two 
successive solutions, and M nuc is the hydrogen nuclear 
burning rate and M w ; n d is the o ptically-thin wind mass- 
loss rate (see Equation (24) in Kato & Hachisu 1994, 
for more detail). The method and numerical technique s 
are essentially the same a s those inlKato et al. I (20111 ). 
We used OPAL opacities (jlglesias fe Rogers 1119961 ). The 
WD radius (the bottom of hydrogen shell-burning) is as- 
sumed to be the Chandrasekhar radius. The mixing- 
leng th parameter of convect ion a is assumed to be 1.5 
(see iKato fe Hachisu I (2009) for the dependence of the 
light curve on a). Internal structures of the envelope are 
essent ially the same as those in Figure 7 of IKato et al. I 
(|2011l ). We calculate optical and UV light curves from 
the blackbody spectrum with the photospheric tempera- 
ture, T p h- To calculate V magnitude, we use the standard 
Johnson V bandpass and add a bolometric correction of 
0.17 mag (see Section [4j. 

In our model, we simply assume uniform chemical com- 
position of the envelope. PU Vul does not show any 
CO/Ne enrichment but the overall chemical composition 
is almos t consistent with being solar ; slightly subsolar 
of iron (IBelvakina et al. 111984 119891) and helium over- 
abun dance (jAndrillat fe Houziaux 111994 ILuna fe Costal 
120051 ) are reported. Thus, we assume four different sets 
of chemical composition (X, Y, Z) by weight for hydro- 
gen, helium, and heavy elements of the envelope as (0.7, 
0.29, 0.01), (0.7, 0.28, 0.02) (0.5, 0.49, 0.01), and (0.5, 
0.49, 0.006). Here Z = 0.01 is closer to the recent esti- 
mate of hea yy element abun dance of solar composition 
(Z = 0.0128: iGrevesse 1120081 ). The WD mass is assumed 
to be 0.6 M Q as listed i n Table CD Model 4 in Tabled] is 
the same as Model 2 in lKato et al~l (|2011l ). 

A typical classical nova shows heavy element enrich- 
ment (C, O, and Ne) in its ejecta, wh ich is interpreted in 
term s of dr edge- up of WD material (jPrialnik fe Kovetzl 
11984 119951 ). PU Vul shows no indication of such en- 
hancement in spectra, which suggests that the WD is 
not eroded during and before the outburst. The theoret- 
ical model described in lKato et al. I (|2011l) showed that 
only a small part of the accreted matter was lost in the 
optically-thin wind, and the rest was burned to helium 
due to hydrogen nuclear burning and accumulated on the 



WD. Therefore, the WD develops a helium layer under- 
neath the newly accreted material. In the next outburst, 
a part of the helium layer will possibly be dredged up and 
mixed into the upper hydrogen layer. In such a case the 
envelope will become helium-rich like in Models 3 and 4 
in Tableffl 

There are observational evidences of wind mass-loss 
from WDs in som e symbiotic stars. For PU Vul, 
iTomov et al. I (|199lD found broad emission wings in HI, 
He I, He II and N IV lines as well as violet-shifted P 
Cygni type absorption components in H I and He I lines 
in the optical spectra taken in 1990- 91, which they at - 
tributed to the hot component winds. iSion et al. I (|1993l ) 
discussed the onset of Wolf-Rayet type wind outflowing 
from the hot component based on the IUE high resolu- 
tion spectra of 1989-1991, and estimated an upper limit 
of M wi „d < 1CV 5 M Q yr" 1 . For AG Peg, the outburst 
lasted about 150 yr, which suggests a low mass WD 
with no optically-thick winds. The wind mass-loss rate 
from the hot compone nt was estimated to be of th e or- 
der of 10~ 7 M yr- 1 (|Vogel fe Nussbaumer 111991 and 
10 ~ 6 M yr" 1 (jKenvon et aT1ll993| ). The intensity of 
the wind diminished in step with the hot component lu- 
minosity during the decline of the outburst. For AE 
Ara, the wind mass loss rate is estimated to be a few 
times 10~ 8 - lO^ Mp, yr" 1 and the WD m ass to be 
M h sin i ~ 0.4 M Q (|Mikolaiewska et al. II2003T) . 

With such poor information on mass-loss rates, we 
simply assume that an optically-thin wind begins to 
blow when the photospheric temperature rises to log T pn 
(K) ~ 4.0 and the wind continues until log Tph (K) 
~ 5.05 at various rates listed in Table [TJ (e.g., M win d = 
5.0 x 10~ 7 A/ Q yr -1 for Model 1). After the temperature 
reaches log T pn (K) ~ 5.05, the wind mass-loss rate drops 
to M wind ~ 1.0 x 1O- 7 M y^ 1 . 

We cannot accurately determine the WD mass of PU 
Vul o nly from our light curve analysis. IKato et al. I 
(|2011| ) obtained a plausible range of the WD mass, 0.5 
- 0.72 M Q corresponding to a reasonable range of the 
wind mass-loss rates. In the present paper, we adopt an 
0.6 M WD as a standard model of PU Vul (see Section 
12.21 for more detail). 

2.2. Continuum UV Light Curve 

In classical novae, a narrow spectral region around 
1455 A is known to be emission-li ne free and can 
be a representative of continuum flux (|Cassatella et al.1 
2002) . This continuum band has been used to determine 
distances to several class i cal novae (lHachisu fe Kato~1 
I2006t lHachisu et al. I l2l)08t IKato et al. II2009D and also 
used in analysis of PU Vul (|Kato et al. 1120111 ). As PU 
Vul shows much weaker emission lines in its spectra than 
classical novae, we can use three other wavelength bands 
around 1350, 1490 and 1590 A of a 20 A width, in 
addition to the UV 1455 A band. Figure [1] depicts light 
curves of these four narrow bands, extracted from the 
IUE data archival 

During the outburst, the photospheric temperature 
gradually rises and the photospheric radius shrinks while 
the bolometric luminosity is almost constant. Thus, a 
UV light curve has the peak at a certain temperature. 

1 http://sdc.lacff.inta.es/ines/ 
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TABLE 1 

Model Parameters 



Subject 


Model 1 


Model 2 


Model 3 


Model 4 


Units 


X 


0.7 


0.7 


0.5 


0.5 




Y 


0.29 


0.28 


0.49 


0.494 




Z 


0.01 


0.02 


0.01 


0.006 




WD mass 


0.6 


0.6 


0.6 


0.6 


Mq 


M bol a 


-5.44 


-5.36 


-5.63 


-5.68 


mag 




-5.61 


-5.53 


-5.80 


-5.85 


mag 


T a 

-^pcak 


4.6 


4.2 


5.5 


5.7 


10 37 erg s _1 


maximum radius c 


63 


60 


61 


64 


Rq 


initial envelope mass 


4.0 


2.6 


3.4 


4.6 


KT 5 Mq 


H-burning rate d 


1.7 


1.6 


2.9 


3.0 


10~ 7 Mq yr" 1 


assumed wind mass-loss rate (T < 5.05) c 


5.0 


3.0 


2.0 


3.0 


lO" 7 Mq yr" 1 


assumed wind mass-loss rate (T > 5.05) f 


1.0 


1.0 


1.0 


1.0 


10~ 7 Mq yr" 1 



Typical values of the optical fiat peak at log T p h (K) —3.9. 
We adopt M boi - 0.17 mag. 



The radius reached before j 



=3.9. 



r ph (K) 

d Values at logT ph (K) =4.5. 

e Optically-thin wind from logT p h(K) = 4 to 5.05. 
Optically-thin wind from logT p h(K) = 5.05 to the end of hydrogen burning 



Figure [T] also shows theoretical UV light curves that rep- 
resent continuum emission in each wavelength. These 
four light curves show basically a similar behavior, be- 
cause each wavelength is close. In a shorter wavelength 
band, the UV flux reaches maximum slightly later than 
in the other longer bands as indicated by upward arrows. 

The flux at the observed peak is obtained to be -F1350 = 
5.7 x 1CT 13 , F1455 = 5.6 x 1CT 13 , F1490 = 6.7 x lCT 13 , 
and ^1590 = 6.1 x 1CP 13 erg s _1 cm~ 2 A _1 , respectively. 
If the emission can be approximated by blackbody, un- 
absorbed peak fluxes should be larger in a shorter wave- 
length band, while the absorbed fluxes are in the inverse 
order. Comparing these peak fluxes, we see that the 
1455 A band flux is too small, because the peak flux is 
more absorbed by cool winds from the M giant compan- 
ion than in the other three bands of 1350, 1490 and 1590 
A (IShore & Aufdenberg |[l99l . The excess of ^1490 may 
be explained by contamination of emission lines. Con- 
sidering these effects, we use the 1590 A band in the 
following discussion. 

Figure [T] also shows model light curves of the 0.6 Mq 
WD with the chemical composition of X = 0.7, Y = 
0.29, and Z = 0.01. Each band light curve is made from 
blackbody emission of our evolution model. Here, we 
assume four optically-thin wind mass-loss rates of 4-8 
xl0~ 7 Mq yr -1 . For a higher mass-loss rate, the evolu- 
tion is faster and the UV light curve shape is narrower. 
All these light curves more or less agree with the obser- 
vational UV light curve in each wavelength band, and we 
chose the 5 x 10~ 7 Mq yr -1 as having the best agreement 
with these data points. 

For a given chemical composition IKato et al. I (|2011[ ) 
obtained a range of the WD mass that reasonably well 
reproduces the UV light curve for reasonable rates of 
the optically-thin mass-loss. The lowest WD mass is ob- 
tained for a very large wind mass-loss rate of 1 x 10 _6 M Q 
yr -1 , while the highest WD mass is for no wind mass- 
loss. For example, if we fix the chemical composition 
to be X — 0.7 and Z = 0.01, a plausible WD mass is 
between 0.52 and 0.72 Mq, corresponding to the wind 
mass-loss rate of 1 x 10~ 6 Mq yr -1 and no mass-loss, re- 
spectively. These ranges of the WD mass are summarized 



in Tabled] for four specified chemical compositions. This 
table also shows a range of the bolometric luminosities 
at the optical fiat peak. The larger the bolometric lumi- 
nosity, the more massive the WD. Combining these the- 
oretical bolometric luminosities with the observed mag- 
nitudes, we can derive a range of the distance moduli, 
(m — M)v, which are shown in the last column of Table 

m 

3. EXTINCTION AND DISTANCE 

Before deriving physical parameters of the nova, 
we must estimate the extinction and distance. The 
reddening was estimated by various methods, H I 
Balmer line ratios, He II emission line ratios, in- 
terstellar optical/UV absorption features, and com- 
parison between the o bserved optical/near-IR spec- 
tra and some sta ndard s (IBelvakina, et al. Ill982bl 119841: 
Friediung et al.1 Il984t iKenvon I 119861 : iGochermann I 
199H iVogel fe Nussbaumer I Il992t 



Rudy et al. Hl999HLuna fe Costa lT2005D . They are un- 
fortunately scattered in a broad range of E(B — V) = 
0.22 — 0.53. Thus, we have made our own estimates based 
on the theoretical light curves (Section |3.1[) and compar- 
ison between spectral classification and colors (Section 



1986; 

Hoard et al. I 119961: 



3.1. Extinction from Model Light Curves 

From our light curve fittings, we get relations on the 
extinction E(B — V) and the distance d to PU Vul. The 
distance modulus is 

(m - M) v = A v + 5 log (d/1 kpc) + 10, (1) 

where A v = R v x E(B - V) and R v = 3.1. In the op- 
tical maximum phase, 1979-1986, except the eclipse, the 
mean magnitude is obtained to be V = 8.59 ± 0.06 (see 
Tabled]), whereas the absolute bolometric magnitude is 
M hoi = -5.44 from Model 1 (Table[T]). Here, we adopt a 
bolometric correction of BC(U)=0.17 (see Section|4]), as 
a representative value for an extended WD photosphere 
during the A-F spectral phase. Then, we have 

14.20 = 3.1x E{B - V) + 5 log (d/1 kpc) + 10. (2) 
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TABLE 2 

Range of Distance Moduli 



Composition 


WD Mass a 


^bol 


^pcak 


-M v c 


(X, Y, Z) 


(Mq) 


(10 37 er 


5 s" 1 ) 






(0.7, 0.29, 0.01) 


. 0.52 - 0.72 


3.4 - 


6.1 


13.87 


- 14.52 


(0.7, 0.28, 0.02) 


. 0.5-0.67 


3.0 - 


5.0 


13.75 


- 14.30 


(0.5, 0.49, 0.01) 


. 0.5-0.62 


3.0 - 


5.7 


14.03 


- 14.45 


(0.5, 0.494, 0.006) 


. 0.53 - 0.65 


4.5 - 


6.5 


14.19 


- 14.58 



a A range of the WD mass obtained from the UV light curve fitting. 
The lower limit corresponds to the case of a very large mass-loss rate of 
1 x Mq yr — 1 , while the upper limit is the extreme case of no wind 

mass- loss fsee lKato et al. II20TTI) . 
b Values at logT (K)=3.90. 

c We adopt V pca k — 8.59 for the optical flat peak. Theoretical absolute V- 
magnitudes are calculated as My — iWbol — 0.17. 




8000 
JD 2,440,000 



(day) 



Fig. 1. — UV light curves for four narrow bands at, a) 1350 A, 
b) 1455 A, c) 1490 A, and d) 1590 A. Theoretical light curves 
are also shown for an assumed distance of 1 kpc (with right-side 
axis) and no absorption; They are 0.6 Mq WDs with the chemi- 
cal composition of the envelope X = 0.7 and Z = 0.01 with four 
different optically-thin wind mass-loss rates; Dash-dotted curves: 
4 X 10 _7 M o yr" 1 . Solid curves: 5 X lO^Mgyr -1 (Model 1). Dot- 
ted curves: 6 X 10 -7 MQvr -1 . Dashed curves: 8 X lO _7 M0yr _1 . 
Short vertical lines in panels b), c), and d) show the period of the 
second eclipse at the optical V band between JD 2,449,270 and 
2,449,610. Among the observational data, the open circles denote 
the ones with low accuracy because they were observed with a short 
exposure time (< 1000 s) or it is obtained from very noisy spectra. 
The red arrows indicate the epoch of the UV maximum in each 
wavelength band. 



This equation gives a relation between E(B — V) and 
d for a specified model, Model 1, which is depicted in 
Figure [2] There is another possible way of fitting. In 
1979, PU Vul showed a spectral type of FO I without 
emission lines, and its magnitude was about V — 8.87 
(see Figure [8]). If we take a bolometric correction typica l 
for F0 I/II, BC(\0=0.13 (jStraizvs fc Kuriliene I Il98l . 
we have a larger distance modulus (m — M)y = 14.44 
for the same Model 1. This case is also plotted in Figure 

m 

We have another distance-reddening relation from the 
UV 1590 A light curve fitting, i.e., 



2.5 logFi59o(obs) : 



-2.5 log F1590 (model) 

-A x + 5 log (d/1 kpc). (3) 



Here A\ = Ri^p E(B — V ) an d we ad o pt #1 590=8. 3 
(jFitzpatrick fc Massa I 12007ft . iSeaton I (|1979l )'s for- 
mula gives a similar value of 7.9. Figure QJd) shows 
-Pi59o(obs) = 6.1 x 10 -13 erg s _1 cm ~ 2 A -1 at the 
UV 1590 A peak, whereas Fi 590 (model) = 1.77 x 10~ 10 
erg s _1 cm ~ 2 A -1 with an assumed distance of d = 1 
kpc. Substituting these values into equation ((3]), we get 
a relation 

6.15 = 8.3 x E(B - V) + 5 log (d/1 kpc) (4) 

for Model 1. Figure [5] also shows Equation Q with two 
additional lines in the both sides which represent a pos- 
sible 15% error in the light curve fitting. This error is a 
summation of the accuracy of the absolute flux calibra- 
tion of IUE (~ 5%) and possible contamination of emis- 
sion/absorption line contribution in the region of 1590 
A, which we assumed to be 10 %. 

Combined these two fittings, i.e., Equations ([2]) and 
g]), we obtain E(B - V) = 0.37 and d = 4.1 kpc, which 
are plotted by a black X-mark (the middle one among 
the three Xs). If we assume a different WD mass, we 
get a different relation between E(B — V) and d, because 
My is different for a different WD mass model. The two 
X-marks in the left /right sides in Figure [2] indicate the 
intersection of the two extreme cases of A/wd =0.52 and 
0.72 Mq, corresponding to the lowest and highest WD 
masses (see Table [2]). 

For different sets of chemical composition, we also 
get different intersections which are shown by different 
symbols in Figure [2j From these points we see that 
E(B - V) = 0.37 is almost independent of the WD 
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0.6 



0.5 - 



0.4 - 



0.3 - 



0.2 - 




1.5 



- 0.5 



0.1 




Distance (kpc) 

Fig. 2. — Distance-reddening relation of PU Vul. The distance-reddening relations obtained from the optical light curve fitting of Model 
1; black solid line: (m — M)y = 14.20, i.e., Equation (0; black dotted line: (m — M)y = 14.44. The red solid line represents the UV 
1590 A light curve fitting of Model 1 (Equation (0). The two red dotted lines beside the line represent possible ±15% errors in absolute 
UV flux. The distance-reddening relation derived from the K-magnitude fitting (Equation JT)) is depicted by the green solid line with ±0.2 
mag error lines in both sides. The central black cross indicates the intersection of the UV and optical fluxes of Model 1. The crosses in 
both sides represent the intersections for J%d = 0.52 (left) and 0.72 Mq (right) WD models with the chemical composition of X = 0.7 
and Z = 0.01, i.e., the lowest and highest WD masses in Table [2] Other models are depicted by the different symbols: blue filled circles for 
Model 2, green open squares for Model 3, red triangles for Model 4. The intersections for the same models but for different extinction law 
are also plotted by blue symbols: (R v , R 1590 )=(2A0, 6.29) for HD 185418, (2.48, 8.03) for HD197702, and (2.91, 8.84) for BD ±35 4258. 
The same symbol indicates the same chemical composition and the middle point of the three same symbols corresponds to the models in 
Table [T] The name of each star is indicated besi de the corresponding group of points with first few letters. The open star mark represent 
d=4.7 kpc and E(B - V) = 0.3. See Section |3~T1 for more detail. 

mass or chemical composition. This is because we use 
the same response (passband) functions to derive My 
and log ^1590 (model) from blackbody spectrum of each 
model, and therefore, the ratios of the two values are 
common in all the models. As a result, these two equa- 
tions yield a common value of E(B — V) independent of 
the model. On the other hand, the distance depends on 
the WD mass and chemical composition (X), because a 
more massive WD/smaller X has a larger photospheric 
luminosity, which results in a larger distance. In this 
way, we could not determine the distance only from the 
light curve fittings. We can constrain the distance corre- 
sponding to a permitted range of the WD mass as listed 
in Tabled! 

It should be noted that Galactic interstellar absorption 
has very large un certainty around the averag e value we 
adopted (see e.g. iFitzpatrick fc Massa 1 12007). Unfortu- 
nately, the extinction curve are not known in PU Vul, 
nor in the stars close to its same sight line. The closest 
stars are relatively far away; towards HD 185418 and HD 
197702 , both ~ 11 deg away from PU Vul, and towards 
BD ±35 4258, ~ 15 deg away. The values of (R v and 
-R1590) are (2.40, 6.29) for HD 185418, (2.48, 8.03) for 



HD 197702, and (2.91, 8.84) for BD ±35 4258. Using 
these values, we obtain the intersection from Equations 
(QJ and ([3]) , which are also shown in Figure [2j These 
extinction estimates strongly depend on the adopted ex- 
tinction curve. 

Thus, we could not accurately determine the distance 
and E(B — V) from the light curve fittings of UV/optical 
bands. 

3.2. Optical Color Excess 

A direct estimate of the color excess E(B — V) can be 
derived from the comparison of observed B — V colors 
and spectral types of the hot component during 1979- 
1983, i.e., at the optical maximum when strong emis- 
sion lines were still absent. In particular, we have calcu- 
lated one-month averaged B — V and U — B colors from 
the p u blished photomet r y dBelvakina et al. Hl982aL 11 985. 
| 1990t iKolotiloy l Il983t iMarerave I 119791: iMahra et al. I 
n979tlWhitnev Hl979l:IBruch 1119801: IC hochol et al. ( ll98lT 
iPurgathofer fc Schnell 1 119*821 ) for several ep ochs when 
the spectral classification of PU Vul was made (iKolotilov 



1983; Bcl vakina et al. I 119841: iMurset fc Nussbaume" 



1994 and references therein). Table [3] shows the one- 
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TABLE 3 
Color excess in optical 



Date 


Sp. Type a 


< V > 


< B-V > b 


(B- V) C 


E(B - V) 


< U-B > b 


(U - B) C 


E(U - B) 


Apr 1979 


A7II 


8.84 


0.39 ± 0.01 


0.12 


0.27 


0.24 ± 0.01 


0.02 


0.25 


May 1979 


FOI 


8.93 


0.42 ± 0.01 


0.19 


0.23 


0.29 ± 0.01 


0.16 


0.13 


May 1979 d 


FOI 


8.98 


0.51 ±0.02 


0.19 


0.32 


0.43 ±0.03 


0.16 


0.27 


Jul 1979 


FOI 


8.84 


0.44 ± 0.01 


0.19 


0.25 


0.34 ±0.01 


0.16 


0.18 


Sep 1979 


F5I 


8.78 


0.57 ±0.01 


0.40 


0.17 


0.37 ±0.02 


0.32 


0.05 


Dec 1979 


F3I 


8.78 


0.69 ±0.01 


0.31 


0.38 


0.50 


0.27 


0.23 


Aug 1981 


F5I 


8.55 


0.75 ±0.01 


0.40 


0.35 


0.47 ±0.01 


0.32 


0.15 


Aug 1981 d 


F5I 


8.64 


0.71 ±0.01 


0.40 


0.31 


0.55 ±0.01 


0.32 


0.23 


Jun 1982 


F5I 


8.50 


0.78 ±0.01 


0.40 


0.38 


0.44 ±0.02 


0.32 


0.12 


Sep 1982 


F5-8I 


8.67 


0.74 ±0.01 


0.49 


0.25 


0.42 ± 0.02 


0.36 


0.06 


Nov 1982 


F3-4I 


8.40 


0.57 ±0.01 


0.33 


0.24 


0.19 ±0.01 


0.28 


-0.09 


Dec 1982 


FOI 


8.43 


0.58 ±0.01 


0.19 


0.39 


0.23 ±0.03 


0.16 


0.07 


Jul 1983 


FOI 


8.39 


0.61 ± 0.01 


0.19 


0.42 


0.27 ±0.01 


0.16 


0.11 



a These spectral types arc taken from Kolotilov (1983). Belyakina ct al. (1984), and Miirsct Sz Nussbaumcr (1994, and 
references therein). 

The average B — V and U — B color indices are calculated using the d ata fro m Belyakin a et al. I I1982al . IT985I . 11990) ; 
IKolotilov I i 19831 ): [Margrave I lH979f ): [Mahra et al. I 1119791 1: Iwhitnev I lH979f ): [Br"uch I Ijl980l). 

c The intrinsic color indexes of supergiants, (B — V)o and (U — B)o, arc from [Straizvs (1992). 

d The UBV colors in May 1979 JChochol et al. H198ll) and in August 1981 dPurgathofer fe Schnell 1119821) show some sys- 
tematic offset with respect to those from other sources, and arc displayed separately. 



month averaged B — V and U — B colors and spec- 
tral classification corresponding to each epoch. Assum- 
ing that the nova envelope of PU Vul had a typical 
spectral type, we can estimate the extinction with the 
comparison to the i ntrinsic color in dex corresponding 
to its spectral type (jStraizvs I fl992f ) . The B-V and 
U — B colors of PU Vul are in good agreement with 
those of bright supergiants for 1979-81, and the av- 
erage ratio of E{U - B)/E{B - V) = 0.66^;^ esti- 
mated for this period agrees with th e predicted value 
of - 0.63 - 0.69 ijCiaxdelli et al~lll989h . Since 1982 there 
is some contamination by the nebular continuum emis- 
sion, especially in the U — B color. We thus obtain 
E{B -V) = 0.30 ± 0.02 for the whole 1979-1983 period, 
and E{B — V) = 0.29 ± 0.02 for 1979-81. 

This value is in good agreement with the color ex- 
cesses estimated from various emission line ratios. We 
estimate E(B — V) ~ 0.2-0.3, from our measurements 
of the He II 1640 emission line fluxes from IUE spec- 
tra obtained in October 1991 and August 1992 com- 
bined with th e optical He II 4686 line flu xes for the 
same epochs (jAndrillat fe Houziaux I 119951 ). assuming 
electron temp erature between 10 000-20 000 K. Similarly, 
IRudv et al. I (p9l derived E(B-V) = 0.22 ± 0.10 f rom 
the O I line ratios whereas iLuna fc Costa I (|2005l ) re- 
ported E(B — V) = 0.29 resulting from the Balmer H I 
line ratios. 

Our extinction estimates agree with estimates for the 
total Galactic extinction towards PU Vul. We esti- 
mated E(B — V) = 0.35 from the Galactic extinc- 
tion distribution based on COBE and IRAS maps com- 
bined with the Leiden -Dwingeloo maps of H I emissions 
(jSchleeel et al. |[l99l . Also E(B - V) ~ 0.2 is obtained 
from the dust map based on analysis of 2MASS photom- 
etry ([Marshall et al. Il2006f ). 

3.3. Distance Estimates from Pulsating RG 

An independent way of distance estimate comes from 
the pulsating RG companion of PU Vul. There is a 
well-known relation between the pulsation period and its 
luminosity for Mira variables, applicable also to semi- 



regular variables pulsating in the fundamental mode, 

M K = -3.51 x (logP(day) - 2.38) - 7.25, (5) 

with an error of ~ 0.2 mag (jWhitelock et a l. 2008). For 
the 218 day pulsation, we get the absolute K magni- 
tude of M K = -7.10 ± 0.2. During the 1980 eclipse 
the a verage K mag change s from 6.26 mag to 6.51 
mag dBelvaki na et al. 1 11985[ ) and its average value is 
< K >= 6.37 ±0.04 mag. 

On the other hand, the difference between the absolute 
and apparent K-magnitudes is written as 

K - M K = 0.353 xE(B-V)+ 5 log (d/lkpc) + 10, (6) 

her e we adopt the reddeni ng law of Ak — 0.353 x E(B — 
V) (jCiardelli et al. 1119891 ). Therefore, we get 

13.47 = 0.353 x E(B - V) + 5 log [d/l kpc) + 10. (7) 

This relation is plotted in FigureO For a particular value 
of E(B - V) = 0.3, we get d = 4.7 kpc. Figure [2] shows 
that the intersections of the three solid lines, i.e., UV, 
optical, and IR, are relatively close, and E[B — V) = 0.3 
from optical color excess is also close to these points. 
This indicates that our various methods are consistent 
with each other. Considering ambiguity of each method, 
we use E(B — V) = 0.3 and d = 4.7 kpc as a reasonable 
estimate. 

4. EVOLUTION OF TEMPERATURE AND RADIUS 

We have deduced the temperature and total luminosity 
of the hot component (WD) of PU Vul using IUE spectra 
as well as published photometry and other useful infor- 
mation. The radius is calculated from the temperature 
and luminosity using the Stefan-Boltzmann law. Here 
we assume E(B - V) = 0.30 and d = 4.7 kp c We 
use the extinction law of iFitzpatrick fc Massa I (|2007l ) 
for the narrow band c ontinuum and emission lines, and 
ICiardelli et al. I (|1989f ) for the broad band optical and 
near infrared photometry, respectively. The resulting 
values are listed in Table [4] For several epochs we have 
been able to calculate them using more than one method, 
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TABLE 4 

Temperature and Radius of Hot Component 3 - 



Date 


JD 


Spectral 


T h 

Ion D 


V 




M ho i 








Method c 




2,400,000+ 


type 






(K) 




(£©) 


(R@) 






04/1979 


43980 


A7 II 




8.84 


7900 


-5.36 


11070 


56 


1 




05/1979 


43991 


F0 I 




8.93 


7400 


-5.23 


9820 


60 


1 




07/1979 


44070 


F0 I 




8.84 


7400 


-5.32 


10670 


63 


1 




12/1979 


44222 


F3 I 




8.78 


6900 


-5.41 


11590 


75 


1 




08/1981 


44834 


F5 I 




8.55 


6500 


-5.66 


14590 


95 


1 




09/1981 


44873 


F5 I 




8.50 


6500 


-5.65 


14510 


95 


2 




06/1982 


45147 


F5 I 




8.50 


6500 


-5.71 


15280 


97 


1 




09/1982 


45229 


F5-8 I 




8.67 


6300 


-5.57 


13430 


97 


1 




11/1982 


45290 


F3-4 I 




8.40 


6800 


-5.79 


16440 


92 


1 




12/1982 


45320 


F0 I 




8.43 


7400 


-5.73 


15560 


76 


1 




07/1983 


45533 


F0 I 




8.39 


7400 


-5.77 


16140 


77 


1 




10/1984 


45989 


A3 I 




8.55 


8900 


-5.77 


16140 


53 


1 




06/1985 


46232 


A3 I 




8.52 


8900 


-5.80 


16600 


54 


1 




09/1986 


46690 


A2 I 




8.72 


9200 


-5.67 


14720 


48 


1 




01/1988 


47176 


_ 


S+ 


_ 


10000 d 


_ 


_ 








06/1988 


47328 




N+ 




15000 d 












10/1988 


47438 




N+2 




35000 


-5.43 


11800 


1. 


3] 


10/1988 


47438 




N+2 




29000 


-5.49 


10704 


4.4 


3a] 


05/1989 


47666 




C+3,N+3 




40000 


-5.47 


12230 


2.3 


3] 


05/1989 


47666 




C+3,N+3 




48000 


-5.47 


15440 


1.8 


3a] 


05-09/1989 


47730 




C+3,N+3 




48000 


-5.70 


15140 


1.8 


1 




04/1990 


47795 




0+3 




55000 


-5.81 


16750 


1.4 


4 




07/1990 


48088 




Ne+2 




41000 d 


-5.64 


14320 


2.4 


4 




1 1 /1 990 


48217 




He+2 




65000 


-5.65 


14500 


0.95 


5 


, SWP40155 


04/1991 


48352 




He+2 




65000 


-5.64 


14350 


0.94 


5 


, SWP41299 


09/1991 


48522 




Hc+2 




67000 


-5.59 


13710 


0.87 


5 


, SWP42536 


10/1991 


48559 




He+2 




70000 


-5.55 


13230 


0.78 


5 


, SWP42937/8 


08/1992 


48858 




He+2 




77000 


-5.67 


14660 


0.68 


5 


, SWP45415 


06/1995 


49886 




Ne+4 




97000 


-5.58 


13540 


0.41 


6 




06/1995 


49886 




Ne+4 




97000 


-5.68 


14860 


0.43 


4 




06/1996 


50237 




He+2 




90000 


-5.31 


10570 


0.42 


5 


, SWP57322/3 


06/1996 


50237 








90000 


-5.42 


11700 


0.36 


4 




09/1996 


50342 




He+2 




83000 


-5.05 


8300 


0.44 


5 


, SWP58251/2 


09/2001 


52154 




Fe+6 




100000 


-4.68 


5940 


0.26 


6 




09/2001 


52154 




Fe+6 




100000 


-5.42 


11700 


0.36 


4 




26/3/2003 


52818 




Fe+6 




165000 


-4.97 


7730 


0.11 


4 




9/4/2004 


53105 




Fe+6 




> 99000 












04/2004 


53110 








150000° 


-4.97 


7720 


0.13 


■1] 


07/2006 


53930 




0+5 




> 114000 












07/2006 


53930 




Q+5 




150000° 


-4.70 


6030 


0.11 


'1] 


06/2011 


55740 








150000° 


-5.02 


8090 


0.13 


'1] 



a We assume d = 4.7 kpc and E(B - V) = 0.30. 
k Highest ionization stage. 

c Methods used in deriving the results (see Scctionf^for details). [1] Supcrgiant phase method; [2] Integration of SED; 
[3] Black body fit to a short wavelength IUE spectrum with T as a free parameter; and [3a] with T from the highest 
ionization stage observed; [4] Miirsct &: Nussbaumcr (1994) method based on UB observations of the nebular phase; 
[5] analysis of He II1640 emission line and ultraviolet continuum; [6] based on He II 4686 emission line fiux. 
d Ionization stage and taken from Miirsct Sz Nussbaumcr (1994). 
e Arbitrary assumed. 



and the differences in the results give an idea about the 
uncertainties of these methods. 

The third column of Table |H shows t h e spectral classifi- 
cation adop ted fromlKolotilov I (119831). IBelvakina et al. I 
(1984), and lMiirset fc Nussbaumer I (|1994 Table 3). and 
the fifth column shows the average of observed V mag- 
nitud es (|Kolotilov I ["19831; IBelvakina et al. iri982al 119851 
1990) of PU Vul during the optical maximum, 1979- 
1986. The outbursting component of PU Vul showed 
spectra similar to that of an early F-type supergiant, 
gradually evolving towards an A supergiant. Since there 
is no strong nebular contribution during this phase, we 
have assumed that the spectral type is an appropri- 
ate indicator for the effective temperature and bolo- 
metric correction , and that the observed V magnitude 
(|Kolotilov I [1981 IBelvakina et~al~l H982aL [19851 [19901) 
represents that of the hot component. We have adopted 



the temperatur es and bolometric correction s for A-F su- 
pergiants from IStraizvs &: Kuriliene I (|1981f) , and adopt 
M0(bol)=4.75 for the absolute bolometric luminosity of 
the Sun. 

The hot component luminosity in September 1981 (JD 
2 444873) have been estimated directly by integrating 
the spectral energy distribution (SED) from ultravio- 
let to infrared. To get the SED we have combined 
the very long exposure IUE spectra (SWP 15110, IWR 
11627 and LWR 11628) from 27 September 1981 with 
IBelvakina et al. 1 (|1985D spectrophotometry performed 
on 25 September 1981 and JHK photometry obtained 
on 22 September 1981. The SED has been corrected 
for the reddening. The resultant bolometric magnitude 
is rnboi.o = 7.71, which derives the absolute bolometric 
magnitude M^ \ — —5.65 with the distance d = 4.7 kpc. 
This value, —5.65, shows an excellent agreement with the 
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potential (IP) observed in the 



( JD 2447438-795) and published op- 
dMunari fc Zwitter I [2002T lYool [20071: 



5000 



5x10^ 



JD 2,440,000+ (day) 



Fig. 3. — Development of the temperature (open triangles and the 
left ordinates) and radius (filled circles and the right ordinates) of 
the hot component (WD) of PU Vul, taken from Table[4] Two red 
upward arrows indicate lower limits of temperature. Model lines 
denote the photospheric temperatures and radii of the 0.6 A/q WD 
with X = 0.7 and Z = 0.01 for five different wind mass-loss rates. 
Dash-dotted lines: 4 X 10~ 7 Mq yr^ 1 . Thick solid lines: 5 X 
10~ 7 Mq yr" 1 (Model 1). Dotted lines: 6 X 10~ 7 M Q yr^ 1 . Thin 
solid lines: 7 X 10~ 7 M yr" 1 . Dashed lines: 8 X 10~ 7 Mq yr _1 . 
Three black downward arrows indicate the central times of the 
1980, 1994, and 2007 eclipses. 

average value, —5.66 in August 1981, derived from the 
observed V mag and spectral type (see Table [4]). 

We have also used this SED to estimate the bolomet- 
ric correction corresponding to this particular date, in 
September 1981. We obtain V a = 7.57 from the SED. 
Thus, the bolometric correction is calcu lated as BC(V) = 
mboi.n — Vo = 0.14. On the other hand, iBelvakina et al. I 
(1985) obtained V = 8.45 on JD 2444873, which is 
corrected to be Vo = 7.52 with an extinction of 3.1 x 
E(B — U)=0.93. Combining this with mboi,o = 7.71, we 
get BC(V) = 0.19. These BC(V) values are somewhat 
larger than BC(V) = 0.08 corresponding to F5 I spectral 
type estimated at this epoch. The difference may reflect 
the lower density in the nova envelope than that in the 
brightest F-type supergiants. In the present work, we 
use BC(V) = 0.17, the mean value of 0.14 and 0.19. 

For the nebular phase, whenever possible, the temper- 
ature of the hot component (WD), Tj,, has been esti- 
mated from the equivalent width of the He II 1640 emis- 
sion line measured from IUE spectra (identified in the 
last column of Tabled]). Although the high resolution 
HST/GHRS spectrum taken in October 1994 (the 1993/4 
eclipse egress) suggest ~ 20% of O I] 1641 line contribu- 
tion to the He II fluxes derived from lower resolution 
spectra, the O I] 1641 line is not visible in the well ex- 
posed high resolution spectra SWP45417 and SWP57730 
taken before and after the eclipse, respectively. There- 
fore, we assume a negligible contribution of O I] to our 
measurements of the He II 1640 line. The luminosity 
of the hot component, Lh, has been calculated from the 
He II 1640 flux, assuming blackbody (Th) and that the 
He II lines are produced by photoionization followed by 
recombination (case B). Th has been also estimated from 
the IUE flux at 1350 A, assuming that it is emitted by 
blackbody (Th). At most epochs in Table2]these two val- 
ues agree with each other, and a mean value is adopted 
for the final Th- 

At several epochs, Th has been derived from the 



highest ionization 
IUE spectra 
tical spectra 

iTatarnikova fc Tatarnikov 201 W. We used the relation 
Th/lO OO ~ IP (eV) found bv iMiirset fc Nussbaumer I 
(|1994ri . For two epochs of JD 2452154 and 2452818, 
Th is derived from the ratios of H/3, He II 4686 
and He I 5876 emission line fluxes published by 
ITatarnikova fc Tatarnikovl (|201lD and iLuna fe Costa I 
(|2005[ ). respectively. 

At two epochs of 1988 and 1989, the hot component 
parameters have been derived by fitting a blackbody to 
the short wavelength (A < 1590 A) part of the spectrum 
obtained by combining SWP34405, SWP34406, and 
SWP34407 for JD 2447438, and SWP36301, SWP36302 
and SWP36304 for JD 2447666. 

After 1996, in the absence of IUE spe ctra, Th has been 
derive d from t he U B mag observed by IBelvakina et al. 1 
(|2000fl and by [Shugarov et al. I (l201ll) w ith th e method 
proposed by IMiirset fe Nussbaumer I (|1994h . This 
method assumes that after subtraction of the contri- 
bution from the RG, the optical magnitudes contain a 
direct contribution from the hot star and an indirect 
contribution from the nebulae. Thus, an accurate es- 
timate of the RG contribution is especially important. 
The RG companion is classified to be a spectral type 
of M6 (Section l5.1|) . so its contribution to UB magni- 
tudes c an be neglected. In fa ct, the B — V and U — B 
colors (jShugarov et al. 1 120111 ) suggest that the contin- 
uum is still dominated by the nebular emission, in agree- 
ment with the optical spectra showing o nly faint fla t 
continuum and strong emission lines (e.g. lYool 12007). 
As described later (Sections [5] and [6] and Figure [8]), 
the summation of the WD and nebular contributions 
dominates the V light curve, although it shows a clear 
~ 0.5 mag pulsation owi ng to the RG. Therefore, we ca n 
safely use the method o f IMiirset fe Nussbaumer I (|1994[ ). 
IMiirset fe Nussbaumer I (| 19941 ) also provided the bolo- 
metric corrections to UBV mag of hot component for a 
wide range of hot component temperature, Th. These 
bolometric corrections were derived by model calcula- 
tions with hot component temperature and nebular den- 
sity as free parameters. Although the RG in PU Vul is 
similar to a Mira component of D-type symbiotic sys- 
tems, we have decided to use the bolometric corrections 
for S-types because the U — B color of PU Vul during 
the nebular phase is similar to that predicted by the 
IMiirset fe Nussbaumer I (|1994[ ) model for S-types, and 
the electron density derived for the nebular phase is sim- 
ilar to the values char acterizing the other S-type systems 
(jLuna fc Costall2005h . 

At JD 2 449 886 and JD 2 452 154, the luminosity of 
the hot component, Lh h as been calculated from the 
published He II 4686 fluxes (lAndrillat fc Houziaux~l ll995l; 
ITatarnikova fc Tatarnikov 11201 ll ). assuming a blackbody 
spectrum with Th and that the He II lines are produced 
by photoionization followed by recombination (case B). 

Our estimated temperature and radius are plotted in 
Figure [3J This figure also shows theoretical models for 
the 0.6 Mq WD with the composition of X = 0.7 and 
Z = 0.01, i.e., the same models as in Figure Q] The 
higher the mass-loss rate, the faster the evolution. All of 
these models show more or less good agreement with our 



PU Vul-evolution of WD, RG, and nebula 



10 yr 



9 yr 



86 82 




38 



37 



36 



5.5 



4.5 



log Lff (K) 



Fig. 4. — Evolution of the hot component of PU Vul in the HR 
diagram. Red open circles are our observational estimates, taken 
from in Table [4] Pairs connected by a line segment indicate the es- 
timates obtained for the same day data but with different methods. 
Observational years (two digits except 1979) are attached beside 
the point. Various types of lines denote loci of theoretical models. 
Solid lines indicate, from upper to lower, Model 4 (red), Model 3 
(green), Model 1 (black), and Model 2 (brown). A dashed line de- 
notes a 0.7 M Q WD with X = 0.7 and Z = 0.01. Four dotted lines 
denote 0.5 Mq WD models with different radius and chemical com- 
positions. From upper to lower, a cold WD (the Chandrasekhar 
radius) with X = 0.5 and Z = 0.02, a cold WD with X = 0.7 and 
Z = 0.01 (upper) which is almost overlapped with a cold WD with 
X = 0.7 and Z = 0.02 (lower). The lowest line denotes a hot WD 
with X = 0.7 and Z = 0.02. Dots represent the epoch when nu- 
clear burning extinguishes. Arrows indicate evolution timescales of 
Model 1; from the beginning to logT p h (K)=4.0, and from logTph 
(K)=4.0 to 5.05, i.e., 8.7 yr and 9.6 yr, respectively. 

observational estimates. It should be noticed that the 
theoretical values are those of a blackbody photosphere. 
Even though, they show good agreement with the IUE 
flux (Figure[T]) and also with the estimates obtained with 
quite different methods (Table 2J. 

Figure S] shows the HR diagram of the hot component. 
This figure shows theoretical tracks of the 0.5, 0.6 and 
0.7 Mq WDs for various chemical compositions. The 
observational estimates are also in good agreement with 
our 0.6 Mq models. 

5. LIGHT CURVE MODELS OF ECLIPSES 

Now we present light curve models of the first (1980) 
and second (1994) eclipses. Here we assume spherical 
shapes of the both components, that the inclination an- 
gle of the orbit is i = 90° , that the RG moves in front of 
the hot component (WD) with a constant velocity V m b, 
and that the RG is radially pulsating with a period of 218 
days and its flux changes in a sinusoidal shape around 
the equilibrium magnitude. We also assume that the ra- 
dius of the RG also varies in a fashion of long-period 
Mira variab les (|Thompson et al 1 120021: iWoodruff et al. I 
12004 120081) . and that the radius varies sinusoidally with 
a phase shift of 0.5 to the flux variation, i.e., the radius 
reaches th e minimum at th e max imum brightness as re- 
ported bv IShugarov et aTH (|2011| ). No accretion disk is 
assumed, because there is no observational indication. 

5.1. The First (1980) Eclipse 

We suppose that the 1980 eclipse is total, in which the 
bloated WD is completely occulted by the pulsating RG 
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Fig. 5. — A close up view of the first eclipse. Our model light 
curve is indicated by the red solid line, which is a summation of 
the eclipsed WD photosphere, a constant nebular emission of V = 
14.0, and the RG photosphere with a sinusoidal oscillation around 

V = 13.6 (green dash-dotted line). The mideclipse on JD 2,444,532 
is indicated by a downward arrow. See text for more details. 

companion. The bottom magnitude of V ~ 13 during the 
eclipse seems to be a bit higher than that of a late type 
M-giant, which suggests the presence of a weak emission 
source which was not occulted. Before going to our model 
construction, we need to examine the magnitude of the 
M-giant companion. 

The spectral classification of the RG companion is 
estimated to be M3-M7, but better estimates are ob- 
tained in longer wavelength bands rather than in the 
optical because of c ontam ination by nebular emission. 
iMurset fc Schmid I (|1999D obtained M6-7, using the 
bands in near IR, i.e. A > 8000 A. iBelvakina et al. I 
(1985) derived a similar type, M6.5, during the 1980 
eclipse. This value is uncorrected for the faint nebula 
(V ~ 14), so there may be some fluctuations by ±1 in 
the spectral type. Therefore, we regard M6 as a reason- 
able average spectral type of the M-giant. 

The magnitude of the RG can be es t imate d from its 
.fsf-band magnitude. IBelvakina et al. 1 ()1985D obtained 
< K >= 6.37 during the eclipse, and its reddening cor- 
rected value is < K >o= 6.26 (see Section 13. 3|) . A 
similar value is o btain ed < K >n= 6.17 ± 0.01 from 
IBelvakina et al. 1 ([2000) and iTatarnikova fc Tatarnikov I 
(j2011|) for an average over 1989-2009 As V - K ~ 7.2 
for an M6 III star, (e.g. lStraizvs1ll992D . we get < V > ~ 
13.4. Assuming E(B — V) = 0.3, the visual magni- 
tude of the giant becomes < V >^ 14.3. This value 
is consistent with an averaged pre-outb urst magnitude 
of V = 14.1 ± _ 0.15 jStephgnsod [l97g) and B = 16.5 
(|Liller fc Liller I H979h . This magnitude < V >~ 14.3 
is much darker than the observed mean magnitude of 

V ~ 13 at the bottom of the eclipse, so we need an 
additional source of emission possibly originated from 
optically-thin plasma such as heated RG cool winds. 

We have constructed an eclipse light curve model, 
assuming that the RG mean magnitude, amplitude of 
the pulsation, and brightness of the additional emission 
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TABLE 5 
Eclipse Light Curve Model 



Subject 


1st eclipse 2nd eclipse 


Units 


mideclipse 


. 4,532 


9,447 


JD 2,440,000+ 


total duration (D) 


. 508 


345 


day 


totality (d) 


. 254 


343 


day 


mean RG magnitude 


. 13.6 


13.6 


mag 


amplitude of RG luminosity 


. 75 % 


65 % 




total amplitude of RG in mag a 


. 2.1 


1.7 


mag 


amplitude of RG radius 


. 7 % 


3 % 




Rc/a 


. 0.246 


0.22 




R h /a 


. 0.070 


0.0007 




nebular emission 


. 14.0 


sec Fig. [6] 


mag 



a V(min)-V(max) 

source are free parameters. Figure [5] shows a close-up 
view of the first eclipse and our light curve model. 

A model light curve, that produces a better fit to the 
observed magnitude data, is shown in Figure [3] We ob- 
tain the total duration of the eclipse D — 508 days, the 
totality d = 254 days, and R*/a = (D + d)/(27rP orb ) 
0.244. where a is the separation of the two stars, P or b is 
the orbital period in units of day, and the asterisk denotes 
the specific radius, because it depends on the timing of 
pulsations at the ingress/egress. The equilibrium radius 
of the pulsating RG is R c = 0.247 a. The RG radius 
is smaller than the equilibrium radius at the second and 
third contacts, i.e., 0.93 and 0.96 times the equilibrium 
radius, we obtain R c /a slightly larger than R*/a. 

A bottom magnitude is obtained as a combination of 
the RG equilibrium magnitude and the nebular emission. 
An equilibrium magnitude of the RG darker than V = 
13.8 does not reproduce the wavy bottom shape, even if 
we assume a very large amplitude of the luminosity. A 
combination of the RG equilibrium magnitude of V ~ 
13.6 - 13.8 and a nebular emission of V ~ 14.0 yield 
a better fitting as shown in Figure O Here we assume 
the RG equilibrium magnitude to be V = 13.6 and the 
nebular emission V — 14.0. 

For the oscillation of the RG radius, a good fitting is 
obtained for amplitudes of the radius AR C /R C = 0.06- 
0.08. For a larger amplitude AR C /R C > 0.1, we cannot 
find a good shape of light curves, because a wavy struc- 
ture appears during the ingress and egress. Here, we 
adopt AR C /R C = 0.07. 

Our fitting parameters are summarized in Table [5) 
showing the mideclipse time, i.e., the time when the RG 
center comes just in front of the WD, the total duration 
of eclipse (D), totality (d), apparent magnitude of the 
RG at its equilibrium state, amplitude of the RG lumi- 
nosity in linear scale ALyjLy, corresponding to the to- 
tal amplitude in magnitude (difference between the max- 
imum and minimum magnitudes), amplitude of the RG 
radius, the ratio of the RG radius to the separation a, and 
the ratio of the WD radius to a. Note that the RG mag- 
nitude at equilibrium is not the arithmetic mean of the 
maximum and minimum magnitudes, because we assume 
a sinusoidal variation in the luminosity (linear scale), not 
i n the magnitude (lo garithmic scale). 

iGarnavich I (|1996[ ) supposed th at this eclipse is par- 
tial be cause of a non-flat bottom. iVogel fc Nussbaumer I 
(1992) explained this non-flat bottom shape as a total 
eclipse, but contaminated by two nebular emissions that 
cause the flux excess in the early half and later half, re- 
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Fig. 6. — A close up vi ew of the second eclip se. Observa- 
tional data are taken from Kolotilov ct al. (1995) (crosses) and 
Yoon & Honcvcutt (2000) (open circles). The red solid line indi- 
cates our composite light curve which is a summation of the four 
components, i.e., the pulsating RG (green dash-dotted line), WD 
(magenta solid line), constant "nebula 1" emission of V = 14.0 
(black solid line), and gradually decreasing "nebula 2" emission 
which is eclipsed by 25 % (brown solid line). The RG is pulsating 
around its mean magnitude of V = 13.6. The mideclipse on JD 
2,449,447 is indicated by a downward arrow. 

spectively. However, the wavy bottom shape in Figured] 
is consistent with the RG pulsation. 

5.2. The Second (1994) Eclipse 

PU Vul began to decline in 1987 from the flat max- 
imum and reached V ~ 11.5 just before the second 
eclipse in 1994. The bottom magnitude of the second 
eclipse is V ~ 11.8 (see Figure O, 1.8 mag brighter than 
that of the 1980 eclipse. This eclipse is considered to 
be total, because the continuum UV radiation, which 
has a WD origin, was disappeared during the eclipse 
(INussbaumer fc Vogel lH996t iTatarnikova fc Tatarnikov I 
I2009D . Therefore, the excess flux (Ay ~ 1.8 mag) is a 
c ontribution of hot nebu l ae. 

INussbaumer fc Vogel I (|1996l ) analyzed UV spectra of 
the hot nebulae and found that highly ionized lines dis- 
appeared during the eclipse and recovered after that, 
whereas low-ionized nebular lines were hardly affected. 
This means that the high excitation lines were emitted 
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from a region close to the WD and the low-ionized neb- 
ular lines are emitted from an outer extended region. In 
other words, the nebulae are also partially occulted. 

Thus, there are three sources of emission: the pulsating 
RG, totally occulted WD, and partially occulted nebu- 
lae. For the RG, we assume a similar model as in the 
first eclipse, i.e., the RG is pulsating around the equilib- 
rium magnitude of V = 13.6 but maybe with different 
amplitudes of the luminosity and radius, which are fit- 
ting parameters. For the WD emission, we take Model 1, 
which is shown in Figure |6] (labeled as WD). We assume 
two sources of nebular emission, one is a constant com- 
ponent (labeled "nebula 1") assumed in the first eclipse, 
i.e., V = 14.0. The other is a decreasing component 
which is partially occulted during the second eclipse (la- 
beled "nebula 2"). Its luminosity and decline rate are 
also parameters in order to obtain the best fit. 

Figure [6] shows the resultant light curve. The ampli- 
tude of the luminosity is determined to be 65 % and that 
of the radius is 3 %. For the nebula 2 component, we 
found that a 25 % occultation of the nebula 2 emission 
yields the best fit. We see that our composite light curve 
represents the temporal change of the optical data. 

These fitting parameters are summarized in Table [SJ 
It is difficult to obtain the WD radius from the light 
curve fitting because the ingress and egress of the second 
eclipse is very steep, which indicates the eclipsed object 
is very small. Therefore, we fixed the WD radius to be 
Rh/a — 0.0007, which corresponds to l.Oi?© for a circular 
orbit of a binary consisting of a 1.0M Q RG and a 0.6Af Q 
WD. This assumption has no effects in determining the 
other parameters. 

5.3. M-giant Pulsation and 3rd Eclipse 

Figure [7] shows a periodic modulation of the V mag- 
nitude, which becomes prominent in the later phase of 
the outburst where the WD component becomes dark. 
This modulation is unclear in the flat maximum except 
the first eclipse (Figures [5]) , because the hot component 
is dominant. We regard that this modulation is caused 
by a pulsation of the RG. 

We obtained the pulsation period to be 218 days, as- 
suming that the period is unchanged from the first eclipse 
until 2010. This 218 day period can reproduce well both 
th e first and second ecli pses as shown in Figures [5] and 
61 iChochol et al I (11998ft obtained a 217 day period and 
Shu garov et al. 1 1)201 lh obtained a 217.7 day period. Our 
value is consistent with these periods. 

Figure [7] shows that one of the minima of the RG pul- 
sation accidentally coincides with the time expected for 
the third eclipse in 2007 indicated by an arrow. This 
narrow dip is not the third eclipse of the WD, because 
the duration is too short, and the WD had already be- 
come very dark in the optical band (see Figured]), and an 
occultation of t he WD hardly changes the total bright- 
ness of PU Vul. iShuearov et al. I (|201lD showed that the 
U magnitude is clearly eclipsed in the third eclipse, but 
the V magnitude is not. Our interpretation is consistent 
with theirs. 

5.4. Radii of Cool/Hot Components 

In Sections 15.11 and 15.21 we have obtained R c /a and 
Rh/a for the first and second eclipses. Assuming a circu- 
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Fig. 7.— Light curve of PU Vul for the period JD 2,447,000- 
2,451,000 (upper part) and JD 2,453,500-2 ,457,500 (lower 
part). Observ ational data are t a ken f rom IKolotilov et al.1 
lYoon i 



3) (crosse s) , 



fe Honevcutt I 120001) (small open circles) , 



IKlein et al. I l| 1994ll (squar es) , Kanarnit su et al. I 119911) (open 
stars), and lliiima I 1T989) (middle size open circles). For the 
lower part, observational data are taken from AAVSO (dots) 
and All Sky Automated Survey (ASAS) (crosses). Downward 
arrows indicate the central times of the eclipses, JD 2,449,447 and 
2,454,362. Short vertical lines indicate epochs of the pulsation 
maxima of the M giant assuming a period of 218 days. 



TABLE 6 

Radii of the Cool and Hot Components 3 



RG mass 


a 


flc(lst) b 


fic(2nd) c 




(Mq) 


(Re) 


(*©) 


(Re) 


(Re) 


3.0 


. I860 


459 


413 


131 


2.0 


. 1670 


111 


370 


118 


1.5 


. 1550 


383 


345 


109 


1.0 


. 1420 


350 


315 


100 


0.8 


. 1360 


335 


301 


95.6 


0.6 


. 1290 


318 


286 


90.9 


0.4 


. 1210 


299 


269 


85.5 



a A circular orbit and a 0.6 Mq WD arc assumed, 
k (1st) means the value obtained for the 1st eclipse 
c (2nd) from the 2nd eclipse 



lar orbit, we calculated the binary separation a from Ke- 
pler's third law, a 3 = G(P orb /27r) 2 (MwD + M RG ). The 
resultant radii of the RG and WD are listed in Table [5] 
as well as a for various RG masses. Here we assume a 0.6 
M WD mass. Estimated RG radii (270 - 460 R ) seem 
to be a little bit larger than those of low mass M giant 
stars, which will be discussed in Section [7] (Discussion) . 

For the hot component, we obtain a radius of Rh/a, = 
0.07 for the first eclipse, that corresponds to i? n ~ 85- 
130 Rq as shown in Table [5] For the second eclipse, we 
have fixed the WD radius to be Rh/a — 0.0007, because 
the ingress and egress are too steep to determine the 
radius. Thus, it is not listed in the table. The steep 
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decline/rise suggest Rh/a < 0.001 which corresponds to 
i?h < 1~2 Rq. We can only say that the radius of the 
hot component reduced by a factor of 100 or more. 

In Section S] we have already shown that the WD pho- 
tosphere had shrunk by about two orders of magnitudes 
between the first and second eclipses, from both the ob- 
servational estimates and theoretical models (see Figure 
[3]). The above radius estimates from the eclipses are very 
consistent with these estimates in Figure [3] This is the 
first time that the shrinkage of a nova WD photosphere 
has been measured by eclipse analysis. 

6. COMPOSITE OPTICAL LIGHT CURVE 

Figure [8] shows an observational light curve of PU Vul 
from the beginning of the outburst until 2010. It also de- 
picts our theoretical composite light curve that consists 
of three components, the WD, RG, and nebulae. The 
WD component, in which we use Model 1, well repro- 
duces the observed UV light curve (Figure [5^) as well 
as the optical light curve until 1989. After 1989, PU 
Vul entered a nebular phase and emission-lines dom inate 
the spectra (jlijima Hl989HKanamitsu et al. 111991ft . Our 
WD model does not include line-emission formed outside 
the photosphere, thus the T^-light curve (red solid line) 
decays much faster than the observed one. For the RG 
component, we assume that the equilibrium magnitude 
is constant, V = 13.6, throughout the outburst as we 
did in the first and second eclipses in Section [5] This 
V = 13.6 is indicated by the green horizontal solid line 
in Figure |8p. 

For the nebular emission, we assume two components: 
one is a constant component of V = 14.0 as depicted 
by the dash-three-dotted line (nebula 1) in Figure |BJd. 
We assumed this component uneclipsed at all during the 
first and second eclipses as shown in Figures [5] and [6l 
We suppose this nebula 1 emission originated from the 
RG cool wind, partially ionized by the radiation from 
the hot component. As this emission is faint, it domi- 
nated the total magnitude only in the first eclipse, so we 
have no information on its magnitude whether it changed 
or not. T herefore, we assumed that t his co mponent is 
constant. iTatarnikova fc Tatarnikov I (|2011[ ) found the 
Raman scattered O VI 6830 line in the optical spectra 
taken in mid 2006 and later. This indicates the presence 
of neutral hydrogen, i.e., the RG cool wind. As the WD 
is still hot, a part of the RG cool wind may be ionized. 
So we reasonably suppose this component is still present. 

The other nebula is originated from the WD, the shape 
of which is represented by the blue dashed line (neb- 
ula 2). This component started at the epoch when the 
photospheric temperature of the WD increased to logT 
(K)=4.0 (open square in Fig ure El). This WD-orig i n com - 
ponent was first discussed by Nussba umer et aLl (1988), 
who concluded that the nebular emission is of WD-origin 
because the relative abundances of C, N, and O are close 
to those of classical novae but different from symbiotic 
stars. This component was eclipsed during the second 
eclipse as i n Figure El 

Recently, IShugarov et al. I (|2011ft reported that all of 
the U, B, and V magnitudes are gradually rising after the 
third eclipse while the mean value of the / magnitude is 
almost constant. This indicates that the WD-origin neb- 
ular component is relatively centrally condensed around 
the WD and, at the same time, widely spread out over 



the orbit. Then the width of eclipse by the RG is so 
wide that a whole period of the orbital phase is partially 
eclipsed as shown by the blue dashed line in Figure |H 
This increase in the brightness (U, B, and V) also sug- 
gests that the hydrogen shell-burning on the WD is still 
on-going. 

7. DISCUSSIONS 
7.1. Comparison with Other Works on Eclipses 

IGarnavich I ()1996l ) estimated the relative size of the 
cool component to be R c /a = 0.28 for the first eclipse and 
0.22 for the second eclipse, assuming symmetric shapes of 
the eclipses. Our corresponding values are R*/a = 0.24 
and 0.22, respectively. The difference in the first eclipse 
is explained from the difference in the totality. Assuming 
that the bottom base line at the first eclipse was V = 11.8 
from the second eclipse, Garnavich get a larger totality 
than ours. Thus, R c /a becomes larger than ours. In 
the second eclipse our R c /a is essentially the same as 
Garnavich's, because the radius oscillation of the RG has 
little effects due to small a mplitude (3 %). 

For the hot component, IGarnavich I ([1996ft estimated 
Rh/a = 0.1 for the first eclipse, which is consistent with 
our value of Rh/a — 0.07, considering difficulty of ac- 
curate fitting with the scattered data. For the second 
eclipse Garnavich's value Rh/a — 0.02 is much larger 
than our value of Rh/a — 0.0007. This difference comes 
mainly from the different data sets. Garnavich used the 
AAVSO data that show slower decline/increase at the 
ingress/egress than those in Figure [6l These AAVSO 
data, however, can be also fitted by a more steep light 
curve that yields Rh/a < 0.01. Therefore, in the both 
cases, we can say that the radius of the hot component 
had decreased at least by a factor of ten between the first 
a nd second ec li pses. 

IGarnavich I (|1996ft concluded that the RG radius 
shrunk by 21 % between the first and second eclipses, i.e., 
from R c /a = 0.28 to 0.22. Our values are much smaller, 
10 % (from R c /a = 0.246 to 0.22), but the shrinkage 
of the radius seems to be real because we cannot find a 
parameter set for the same RG radius between the two 
eclipses. We will discuss the shrinkage of the radius in 
the next subsection. 

The orbital period of PU Vul is estimated from the 
mideclipses of the first and second eclipses to be 4915 
days (13.46 yr) (see Table Ell . The orbital perio d 
was obtained as 4918 ± 8 days (|Kolotilov et al. Ill995ft . 
4900 ± 100 davs ([Nussbaumer fc Vogel II1996D. 49 00 ± 9 
(IGarnavich 1 119961) . 4897 days (jShugarov et al~l 12011ft . 
assuming symmetric shapes of the eclipses. Our anal- 
ysis first includes a radius oscillation and the resulted 
non-symmetric shapes of eclipses. However, these effects 
cause only several days off from the symmetry center be- 
cause of small amplitudes of the radial oscillations. This 
is the reason why our new period is close to the previous 
estimates. 

7.2. Comparison with Other Evolution Calculations 

Figure [4] shows the evolution timescale of Model 1, 
18.3 yr from the beginning of the outburst to logT p h 
(K)=5.05. If we do not include the optically-thin wind 
mass-loss, this becomes 46 yr, and the total duration of 
the outburst, from the beginning to the extinguish point 
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Fig. 8. — Optical and UV light curves of PU Vul. (a) UV light curve. Large open circles denote the IUE UV 1590 A band (the same 
as those in Figure [TJi). The solid curve denotes the UV light curve of Model 1. The scale is in units o f 10~ 13 erg s —1 cm " 2 A -1 , (b) 
Optical light curve. Middle size open circles are V magnitudes during the flat phase listed in Table [4] See Kato ct alH 1 120111 ) for the other 
observational data. The red solid and dotted lines indicate optical and bolometric light curves of Model 1. The large open square indicates 
the epoch at logT p i 1 (K)=4.0 when the optically thin wind of 5 X 10 — 7 Mq yr -1 started in our model. The magenta solid line indicates a 
composite light curve of the WD (red solid line), RG at the mean luminosity V = 13.6 (horizontal green solid line), emission of a constant 
component at V = 14.0 ("nebular 1": blue dash-three-dotted line), and a variable component ("nebular 2": blue dashed line). Downward 
arrows indicate the mideclipses. We suppose that the "nebula 2" is eclipsed by the RG companion after 2002 (denoted by 'wide eclipse'). 

of nucle ar burning, is 130 yr. U sing a hydrodynami- 
cal code lPrialnik fc Kovetz I (|1995l ) calculated multicycle 
nova evolution models for various WD masses and ac- 
cretion rates. For a 0.65 Mq WD and a mass accretion 
rate of 1 x 10 _ 'M©yr -1 , no optically thick wind mass- 
loss arose. They obtained the total duration of the nova 
outburst to be £3t>oi=155-176 yr, depending on the core 
temperature, where is the time during which the 
bolometric luminosity drops by 3 mag. For lower accre- 
tion rates (< 10 -8 M Q yr -1 ) strong optically thick winds 
occur which shorten the total duration. Their total du- 
ration is very consistent with our 0.6 Mq WD model, 
considering the different definition of the end point of a 
shell flash; Our definition is for the hydrogen burning ex- 
tinguish point (depicted by the dot in FigureH]), whereas 
Prialnik & Kovetz' is for <3boi time which comes later 
than our extinguish point, thus gives a longer timescale 
than ours. 

Following the referee's suggestion we discuss the work 



bv I Cassisi et al. I (|1998ft and iPiersanti etHI (|1999t 
I200CO who calculated shell flashes on low-mass WDs us- 
ing a spherical symmetric hydrostatic code with the Los 
Alamos opacity. Cassisi et al.'s (1998) models show that 
the envelope extends only down to logT p h (K)=4.5— 4.7 
for a 0.5 Mq WD with m ass accretion rat e s of 2 and 
4 xlO -8 M yr _x . Also in IPiersanti et al. I (|2000[) . the 
temperature decreases down to log T p h (K) =4.1-4.2 only 
in a few exceptional cases. Such high-temperature shell 
flashes may be observed as a UV flash. In other words, 
these calculations do not represent realistic nova out- 
bursts in which the surface temperature drops to log T p h 
(K) < 4.0 at the optical peak. This suggests that their 
numerical code has some difficulties in calculating real- 
istic nova outburst models. 

It should be pointed out that the above three 
works arc obtained with th e Los Alamos opacity , 
not with the OPAL opacity ((Rogers fe Iglesias I [l992t 
llglesias fc Rogers! 1993. 19961), which has been widely 
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used in stellar evolution codes including nova outbursts. 
We are puzzled by the remark in Cassisi et al. that 'the 
Los Alamos opacities are very similar to t he OPAL opac- 
ities' (see the last sentence of Section 4 in I Cassisi et al. 1 
1998). It is well known that the OPAL opacities have 
a strong peak at logT (K) ~ 5.2, while the Los Alamos 
opacities do not (For compariso n with these opacities i n 
a nova envelope, see Figure 15 in lKato fc Hachisu Ill99l . 
This strong peak causes substantial changes in nova out- 
bursts, e.g., acceleration of optically thick winds, and as 
a results, nova evolutions had significantly changed (e.g., 
co mpare hydrod ynamical calculations of nova outbursts 
in lPrialnikl (|1986[) obtained wi th the Los Alamos opacity 
with lPrialnik fc Kovetzl (| 19951) with the OPAL opacity). 
In a less massive WD (< 0.6M©), no optically thick wind 
is accelerated, but internal structures of the envelope are 
significantly different; a density inversion layer appears 
correspo nding to the peak of the OPAL opacity (see Fig- 
ure 7 in lKato et al. I [20111 ). In order to make a reliable 
outburst model of PU Vul, we need to use the OPAL 
op acity, not th e Los Alamos opacity (see also Discussion 
in iKatol [2011) . 

7.3. Pulsating RG Companion 

As described in the previous subsection, our eclipse 
analysis shows that the RG radius decreased by ~ 10 % 
between the first and second eclipses. This radius may 
not be the photospheric radius of the RG defined in near 
IR bands, but the radius of a thick TiO atmosphere, 
which is transparent in K-bemd but opaque in U-band. 
In the pulsating RG atmosphere, the temperature de- 
creases in the expanding phase, which accelerates TiO 
molecule formation, resulting in a large opacity in the 
optical region, which causes a deep minimum in the op- 
tical light curve. The radius, that we obtained from the 
eclipses in U-band, corresponds to the radius of the TiO 
at mosphere. We call this t he "visual photosphere" af- 
ter EeiJ&iGoIdston] ((2002). This radius could be much 
larger than the photospheric radius usually defined with 
-fT-band observation. Therefore, it is very likely that our 
visual photosphere in Table [6] is larger than the RG ra- 
dius in X-band. 

As shown in Section 15.31 the pulsation period of 218 
days had not changed between the first and second 
eclipses. The unchanged pulsation period means that 
the internal structure of the RG had not changed, so the 
if-band photospheric radius should be the same. On the 
other hand, our analysis clarified that the pulsation am- 
plitude in F-band decreased from 75 % to 65 %, and 
also the amplitude of the radius decreased (see Table [5]). 
This suggests that the radius of the visual photosphere 
decreased as the amplitudes of the luminosity and ra- 
dius had decreased. This can be understood as follows; 
The TiO atmosphere is pushed outward in an expanding 
phase, and it pushed far outward when its amplitude is 
larger. Therefore, a larger amplitude results in a larger 
visual photosphere. 

We can estimate the RG radius, using the period- 
luminosity (PL) relations of Mira/semi-regular variables. 
The bolometric luminosity of LMC Mira variables follows 
a PL relation of 

m bo i = (-3.06 ± 0.26) logP + (21.50 ± 0.61), (8) 

where P is the pulsation period in units of day 



(|Glass et al. I 120031) . The fundamental pulsation 
mode of Mira variables corresponds to this sequence. 
With the distance modu lus of LMC, 18.39 ± 0.05 
(jvan Leeuwen et aT~l 120071 ) . we obtain Mboi = —4.05 
for P = 218 day. Therefore, if the RG companion fol- 
lows the PL relations of Miras, its absolute luminosity 
is Mboi = —4.05, i.e., logL = 3300 Lq. Photometric 
studies on a large number of stars indicate another PL 
relation, parallel to the above relation, but about one 
magnitude brighter. This sequence corresponds to the 
first overtone of pulsation in semi-regular variables. In 
this case, we have Mboi = —5.05, i.e., logL = 8300L©. 

The spectral type of the companion is estimated to be 
M6 (see Section 15.11) . The temperature calibration for 
late M giants is relatively well establish ed, and various 
group s give similar values. In particular, iRichichi et al. I 
(1999) give T cS = 3240±75 K an d 3100±80 K for M6 and 
M7 giants, respectively, whereas I Van Belle et al. I (1999) 
report 3375 ± 34 K and 3095 ± 29 K for M6 and M7, re- 
spectively. So, Toff = 3200 ± 100 K for the M6-7 giant in 
PU Vul seems very reasonable. The radius then becomes 
R = 187 ± 12i?0. For the second PL relation we have 
R = 296 ± 19i?0. Comparing these radii with the ones 
in Table O we may conclude that the pulsation of the 
RG companion is consistent with the fundamental mode 
rather than the first overtone, because the visual photo- 
sphere is much larger than the stellar radius (R c ~ 1.8 
times the stellar r adius, that is, R c ~ 330 Rq for the 
fundamental mode: lR"eid fc Goldston ||2002D . 

If the RG pulsation is in the first overtone, the abso- 
lute magnitude is about 1 mag brighter than in the fun- 
damental mode as described above, thus the distance is 
1.6 times larger, i.e., d = 4.7 kpc xl.6 = 7.4 kpc. Such a 
large distance is inconsistent with the optical light curve 
fittings, because E(B — V) becomes too small or nega- 
tive ( see black solid / dotted lines in Figure [5]) , thus we 
cannot construct a consistent model among the optical, 
UV 1590 A, and extinction. Note that the "IR" line in 
Figured! i.e., Equation §5§ is for the fundamental mode 
and the corresponding line for the first overtone is in the 
right outside of the figure. Therefore, we may conclude 
that the pulsation is a fundamental mode. 

Next, we estimate the RG mass using a theoretical re- 
lation obtained from radial-pulsations. It is well known 
that the pulsation constant, Q — P^M/R 3 , is insensi- 
tive to stellar structure, here M is the stellar mass in 
units of Mq, R the radius in units of Rq, and P the pul- 
sation period in units of day. Therefore, the pulsation 
mass is given by 



Numerical calculations show that Q = 0.06 - 0.08 for 
the fundament al mode, and Q = 0.03 - 0.04 for the first 
overtone fe.g.. Kiong. fc Deng 1 12007). Using the radius 
estimated above and P — 218 day, we can estimate the 
RG mass (pulsation mass) to be M = 0.5 - 0.9 Mq for 
both of the fundamental and first overtone modes. The 
0.8 Mq is consistent with the visual photospheric radius 
of 335 Rq in Table [5] for the fundamental mode. 

A different way to estimate the RG mass comes from 
binary evolution theory. A WD of mass ~ 0.6 Mq cor- 
responds to a ~ 2.0 Mq zero-age main-sequence star in 
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the initial-final mass relation derive d from observation 
(e.g., Table 3 in iWeidemannl I2000D . or a ~ 3 Mp, in 
stella r evolution calculation for binaries (jUmeda et al. I 
1999). Then, the companion star should be smaller than 
2 — 3 Mq, because the more massive component in a bi- 
nary evolves first. These values are consistent with the 
above estimate derived from the pulsation theory. 



7.4. X-ray observation 

PU Vul becomes a supersoft X-ray source in the 
later phase of the outburst when the surface tempera- 
ture of the WD b ecomes high enough to emit X-rays. 
iKato et al. I ()2011h estimated the supersoft X-ray flux, 
but it was very uncertain because the long-term evolution 
of the temperature depends on the assumed optically- 
thin mass-loss rate as well as the possible absorption 
due to the RG cool winds. In the present work, we 
confirm that the nuclear burning still continues and 
thus the WD currently evolves toward a supersoft X- 
ray phase. We also showed that the binary system 
contains two different origins of nebulae, i.e., nebula 1 
comes from the RG cool-wind and nebula 2 from the 
WD hot- wind. This cool- wind origin nebula absorbs 
a part of the supersoft X-ray flux, because the neb- 
ula i s partially neutral (Raylei g h sca ttering in 1991- 
1993 lTatarnikova fc Tatarnikov I (120091): Raman scat- 
tered O VI lines since 2006: iTatarnikova fc Tatarnikovl 
(|2011l) ). Therefore, the supersoft X-ray flux should vary 
with the binary phase, i.e., the flux is minimum when 
the RG is in front of the WD and maximum when the 
WD is in front of the RG. The UBV light curves of PU 
Vul (jShugarov et al. |[2011D show such a long term vari- 
ation with the orbital phase. As mentioned in Section [6] 
we explain this variation as an eclipse of nebula 2 by the 
RG (and also possibly by the nebula 1). Therefore, the 
supersoft X-ray flux may also show a similar long-term 
variation. We expect that the X-ray flux will be max- 
imum when the UBV flux is maximum. If the orbit is 
circular, the next maximum will be in June 2014, and it 
is a good chance to detect supersoft X-rays. 

SMC 3 is a symbiotic star consisting of a massive 
WD and an M-giant, which attracts atten tion in relation 
to th e progenitor of type la supernovae (jHachisu et al.l 
2010). Its supersoft X-ray flux and B-magnitude 
show similar long-t erm variations in the orbital phase 
(|Sturm et al. 1 120 111) . Sturm et al. analyzed the X-ray 
variation, but could not give a definite explanation about 
the X-ray variability. We could, however, explain this X- 
ray and B-magnitude variations in terms of wide eclipses 
because it is similar to the UBV variations of PU Vul. 
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8. CONCLUSIONS 

Our main results are summarized as follows: 

1. We present new estimates of the temperature and 
radius of the hot component from a very early phase of 
the outburst (1979) until 2011. These are very consistent 
with our theoretical model of outbursting WDs based 
on thermonuclear runaway events without optically thick 
winds. 

2. We analyzed the first (1980) and second (1994) 
eclipses, assuming sinusoidal variations of the brightness 
and radius of the RG. Both of the eclipses are explained 
as a total eclipse of the WD occulted by the pulsating 
RG. Between the first and second eclipses, both of the 
components shrank in size. The radius of the hot com- 
ponent decreases from ~ 100i?© to ~ O.IRq, which is 
very consistent with our theoretical model. 

3. We are able to construct a composite optical light 
curve that consists of four components of emission, i.e., 
the WD photosphere, hot nebulae surrounding the WD, 
RG photosphere, and nebulae possibly originating from 
the RG cool winds. 

4. We have estimated the extinction and distance with 
various methods, that is, the light curve fittings of optical 
and UV 1590 A bands based on our WD model, direct es- 
timates of color excess, and using i\~-magnitude and P-L 
relation of the pulsating RG companion. These different 
methods yield consistent values of E(B — V) ~ 0.3 — 0.4. 
and d = 4 - 5 kpc. We adopt E(B - V) = 0.3 and 
d = 0.47 in the present paper as representative values. 

5. We interpret the recent long term evolution of V 
magnitude in terms of eclipse of the hot nebula surround- 
ing the WD: the V magnitude gradually decreased from 
2002 and reached a minimum in 2007 and is now in a re- 
covering phase in 2012. This means that hydrogen burn- 
ing is still ongoing. Therefore, we suggest X-ray obser- 
vations around June 2014 to detect supersoft X-rays. 

The authors wish to express great thank to S. Shugarov 
and A. Tatarnikova for providing photometric data as 
well as for discussion on recent observations of PU 
Vul. We are very grateful to M. Takeuti and N. Mat- 
sunaga for providing information and valuable discussion 
of pulsating red giants. We also thank to the anony- 
mous referee for useful comments that helped to im- 
prove the manuscript. We also thank A. Cassatella and 
R. Gonzalez-Riestra for useful discussion, the American 
Association of Variable Star Observers (AAVSO) and 
All Sky Automated Survey (ASAS) for archival data 
of PU Vul. This research has been supported in part 
by the Grant-in-Aid for Scientific Research (20540227, 
22540254) of the Japan Society for the Promotion of Sci- 
ence and by the Polish Research Grant No. N203 395534. 



Belyakina, T.S. et al. 1989, A&A, 223, 119 

Belyakina, T.S. et al. 1990, Bulletin of the Crimean Astrophysical 

Observatory, 81, 28 
Belyakina, T.S., Burnashev, V., I., Gershberg, R.E., Efimov, J.S., 

Shakhovskoi, N.M., & Shenavrin, V.I. 2000, Izv. Krymsk. Astr. 

Obs. 96, 18 

Bensammar, S., Friedjung, M., Chauville, J., & Letourneur, N. 

1991, A&A, 245, 575 
Bruch, A. 1980, IBVS, No. 1805 
Cahn, J.H., 1980, Space Sci. Rev., 27, 457 



16 



Kato et al. 



Cassatella, A., Altamore, A., & Gonzalez-Riestra, R. 2002, A&A, 
384, 1023 

Cassisi, S., Iben, I. Jr., & Tornambe, A. 1998, ApJ, 496, 376 
Chochol, D., Hric, L., & Papousek, J. 1981 IBVS, No. 2059 
Chochol, D., Pribulla, T., & Tamura, S. 1998 IBVS, No. 4571 
Cardelli, J.A., Clayton, G.C., & Mathis, J.S. 1989, ApJ, 345, 245 
Fitzpatrick, E.L., & Massa, D., 2007, ApJ, 663, 320 
Friedjung, M., Ferrari-Toniolo, M., Persi, P., et al. 1984 in The 

Future of Ultraviolet Astronomy Based on Six Years of IUE 

Research, NASA CP-2349, eds. J.M. Mead, R.D. Chapman, & 

Y. Kondo (NASA, Washington,DC) p.305 
Garnavich, P. M. 1996, JAAVSO (The Journal of the American 

Association of Variable Star Observers), 24,81 
Glass, I. S., Evans, T. L. 2003, MNRAS, 343, 67 
Gochermann, J. 1991, A&A, 250, 361 

Gonzalez-Riestra, R., Cassatella, A. & Fernandez-Castro, T. 

1990, A&A, 237, 385 
Grevesse, N. 2008, in Comm. in Astcroseismology, 157, 156 
Gromadzki, M., Mikolajcwska, J., Whitelock, P. Marang, F. 

2009, Acta Astr. 59, 169 
Hachisu, I., & Kato, M. 2006, ApJS, 167, 59 -80. 
Hachisu, I., & Kato, M. 2007, ApJ, 662, 552 
Hachisu, I., & Kato, M. 2010, ApJ, 709, 680 
Hachisu, I., Kato, M., Nomoto, K. 2010, ApJ, 724, L212 
Hachisu, I., Kato, M., & Cassatella, A. 2008, ApJ, 687, 1236 
Hoard, D.W., Wallerstein, G., & Willson, L.A. 1996, PASP, 108, 

81 

Hric, L., Chochol, D. & Grygar, J. 1980, IBVS, No. 1835 
Iijima, T. 1989, A&A, 215,57 

Iglesias, C.A., & Rogers, F.J. 1993, ApJ, 412, 752 
Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 
Kanamitsu, O., Yamashita, Y., Norimoto, Y., Watanabe, E., & 

Yutani, M. 1991, PASJ, 43, 523 
Kato, M. 2012, in Binary Paths to the Explosion of Type la 

Super nova, IAU Symp. 281 eds. R. Di Stefano & M. Orio, in 

press l|arXiv:1110.0055|l 
Kato, M., & Hachisu, I., 1994, ApJ, 437, 802 
Kato, M., & Hachisu, I. 2009, ApJ, 699, 1293 
Kato, M., Hachisu, I., & Cassatella, A. 2009, ApJ, 704, 1676 
Kato, M., Hachisu, I., Cassatella, A., & Gonzalez-Riestra, R. 

2011, ApJ, 727, 72 
Kato, M., Hachisu, I., Kiyota, S., & Saio, H., 2008, ApJ, 684, 1366 
Kenyon, S.J. 1986, AJ, 91, 563 

Kenyon, S.J., Mikolajcwska, J., Mikolajewski, M.,Polidan, R. S., 

Slovak, M. H. 1993, AJ, 106, 1573 
Klein, A., Bruch, A., & Luthardt, R. 1994, A&AS, 104,99 
Kolotilov, E. A., & Belyakina, T.S. 1982, IBVS, No. 2097 
Kolotilov, E. A. 1983, soviet astronomy, 27,432 
Kolotilov, E. A., Munari, U., & Yudin, B.F. 1995, MNRAS, 275, 

185 

Kozai, Y. 1979a, IAU Circ.No. 3344 
Kozai, Y. 1979b, IAU Circ.No. 3348 

van Leeuwen, F., Feast, M. W., Whitelock, P. A., Laney, C. D. 

2007, MNRAS, 379, 723 
LiUer, M., & Liller, W. 1979, AJ, 84, 1357 
Luna, G.J.M., & Costa, R.D.D. 2005, A&A, 435,1087 
Mahra, H.S., Joshi, S.C., Srivastava, J.B., & Dhir, S.L. 1979, 

IBVS, No. 1683 
Margrave, T.E. 1979, IAUC, No. 3421 

Marshall, D. J., Robin, A. C, Reyl C, Schultheis, M., & Picaud, 

S. 2006, A&A, 453, 635 
Mikolajewska, J., Quiroga, C, Brandi, E., et al. 2003, in 

Symbiotic Stars Probing Stellar Evolution, eds, R. L. M. 

Corradi, J. Mikolajewska, & T. J., Mahoney, ASP Conference 

Series, 303 (ASP:San Francisco), p.147 
Munari, U., & Zwitter, T. 2002, A&A, 383, 188 



Miirset, U., & Nussbaumer, H. 1994, A&A, 282,586 
Miirset, U., & Schmid, H.M. 1999, A&AS, 137,473 
Nussbaumer, H., & Vogel, M. 1996, A&A, 307,470 
Nussbaumer, H., Schild, H., Schmid, H.M., & Vogel, M. 1988, 
A&A, 198, 179 

Piersanti, L., Cassisi, S., Iben, I. Jr., & Tornambe A. 1999, ApJ, 
521, L59 

Piersanti, L., Cassisi, S., Iben, I. Jr., & Tornambe A. 2000, ApJ, 

535, 932 
Prialnik, D. 1986, ApJ, 310, 222 
Prialnik, D., & Kovetz, A. 1984 ApJ, 281, 367 
Prialnik, D., & Kovetz, A. 1995, ApJ, 445, 789 
Purgathofer, A., & Schnell, A. 1982, IBVS, No. 2071 
Purgathofer, A., & Schnell, A. 1983, IBVS, 
Reid, M. J., & Goldston, J. E. 2002, ApJ, 568, 931 (Erratum: 

2002 ApJ, 572, 694) 
Richichi, A., Fabbroni, L., Ragland, S., & Scholz, M. 1999, A&A, 

344, 511 

Rogers, F.J., & Iglesias, C.A. 1992, ApJS, 79, 507 

Rudy, R. J., Meier, S. R., Rossano, G. S., Lynch, D. K., Puetter, 

R. C, and Erwin, P. 1999, ApJS, 121, 533 
Savage, B. D., & Mathis, J. S. 1979, ARAA, 17, 73 
Schlegel, D. J., Finkbeiner, D. P. & Davis, M. 1998, ApJ, 500, 525 
Seaton, M. J. 1979, MNRAS, 187, 73 
Shore, S.N., & Aufdenberg, J. P. 1993, ApJ, 416,355 
Shugarov et al. 2011 in "Asiago Workshop on Symbiotic Stars", 

eds. A. Siviero, & U. Munari, Baltic Astronomy, special issue, 

in press 

Sion, E. M., Shore, S.N., Ready, C. J., Scheible, M. P. 1993, AJ, 

106, 2118 
Skopal, A. 2006, A&A, 457, 1003 
Stephenson, C.B. 1979, IAUC 3356 

Straizys, V. 1992, Multicolor stellar photohetry, (Tucson: Pachart 
Pub. House) 

Straizys, V., & Kuriliene, G. 1981, Ap&SS, 80, 353 

Sturm, R., Haberl, F., Greiner, J., Pietsch, W., La Palombara, 

N., Ehle, M., Gilfanov, M., Udalski, A., Mereghetti, S., & 

Filipovi, M. 2011, A&A, 529, 152 
Tatarnikova, A. A. & Tatarnikov, A. M. 2009, Astronomy 

Reports, 53, 1020 
Tatarnikova, A. A. & Tatarnikov, A. M., Esipov, V.F., Tarasova, 

T.N., Shenavrin, V.I., Kolotilov, E.A., & Nadzhip, A.E. 2011, 

Astronomy Reports, 55, 896 
Thompson, R. R., Creech-Eakman, M. J., van Belle, G. T. 2002, 

ApJ, 577, 447 

Tomov,T., Zamanov, R., Iliev, L., Mikolajewski, M., & Georgicv, 

L. 1991, MNRAS, 252,31 
Umeda, H., Nomoto, K., Yamaoka, H., & Wanajo, S. 1999, ApJ, 

513, 861 

Van Belle, G.T. et al. 1999, ApJ, 117, 521 
Vogel, M., & Nussbaumer, H. 1992, A&A, 259, 525 
Vogel, M., & Nussbaumer, H. 1994, A&A, 284, 145 
Weidemann, V. 2000, A&A, 363, 647 
Wenzel, W. 1979, IBVS, No. 1608 

Whitelock, P.A, Feast, M.W., & van Leeuwen, F. 2008, MNRAS, 
386, 313 

Whitney, C.A. 1979, IAUC, No. 3348 
Wood, P.R. 2000, Publ. Astr. Soc. Australia, 17, 18 
Woodruff, H. C, Eberhardt, M., Driebe, T., et al. 2004, A&A, 
421, 703 

Woodruff, H. C, Tuthill, P. G., Monnier, J. D., et al. 2008, ApJ, 
673, 418 

Xiong, D.R., & Deng, L. 2007, MNRAS, 378, 1270 
Yoo, K-H. 2007, Journal of the Korean Astronomical Society, 40, 
39 

Yoon, T.S., & Honeycutt, R.K. 2000, PASP, 112, 335 



